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Abstract 



In this paper we present a simple strategy for the elimination of the translational kinetic 
energy contamination of the total energy in pre-Born-Oppenheimer calculations carried out 
in laboratory-fixed Cartesian coordinates (LFCCs). The simple expressions for the coordi- 
nates and the operators are thus preserved throughout the calculations, while the mathemat- 
ical form and the parametrisation of the basis functions are chosen so that the translational 
and rotational invariances are respected. The basis functions are constructed using explicitly 
, P *. correlated Gaussian functions (ECGs) and the global vector representation. 

First, we observe that it is not possible to parametrise the ECGs so that the system 

> ■ 

, is at rest in LFCCs and at the same time the basis functions are square-integrable with a 

O ■ 

' non-vanishing norm. Then, we work out a practical strategy to circumvent this problem by 

mi 

making use of the properties of the linear transformation between the LFCCs and trans- 
lationally invariant and center-of-mass Cartesian coordinates as well as the transformation 
properties of the corresponding basis function parameter matrices. By exploiting these for- 
mal mathematical relationships we can identify and separate the translational contamination 

rS 

terms in the matrix representation of the kinetic energy operator in the LFCC formalism. 

We present numerical examples for the translational contamination and its elimination 
for the two lowest rotational energy levels of the singlet hydrogen molecule, corresponding 
to para- and ortho-H2, respectively, treated as four-particle quantum systems. 
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1 Introduction 



Expressions for the calculation of intrinsic properties of molecules should be free of 
contributions from the overall translation of the system. In the commonly introduced 
Born-Oppenheimer (BO) approximation, the nuclei are fixed, and thus the transla- 
tional contribution is automatically separated. In several combined electron-nuclear 
orbital approaches [THS] the translational dependence is eliminated automatically by 
fixing one or a few heavy particles. Here, we consider molecules as many-particle 
quantum systems with electrons and nuclei both treated as quantum particles on equal 
footing in the pre-Born-Oppenheimer (pre-BO) quantum theory. 

Traditionally, in rovibrational calculations, in which all nuclei are treated as quan- 
tum particles on a potential energy surface [1HH], the first step is the separation of the 
Cartesian coordinates of the center of mass followed by the definition of a body-fixed 
frame, orientational angles, and internal coordinates. This approach results in the 
replacement of the original laboratory-fixed Cartesian coordinates with curvilinear co- 
ordinates and the corresponding very complicated, translationally invariant rotational- 
vibrational Hamiltonians, see for example [9]. 

Less complicated translationally invariant Hamiltonians are used in full pre-BO 
calculations, where the original laboratory-fixed Cartesian coordinates are replaced by 
some set of translationally invariant Cartesian coordinates (TICCs) and the center-of- 
mass Cartesian coordinates (CMCCs) are separated [TOl - TH] . Although the resulting 
TICCs are rectilinear coordinates, the corresponding Hamiltonian is still complicated. 
It is therefore reasonable to ask whether we can make our calculations even simpler, 
using the original laboratory-fixed-Cartesian-coordinate (LFCC) formalism without 
having to rely on any coordinate transformation at all. In this work we therefore 
explore the usage of LFCCs in pre-BO calculations. 

If LFCCs are used one has to make sure that the energy of the overall translation 
of the system is eliminated. The most straightforward way is the subtraction of the 
kinetic energy operator of the center of mass from the total Hamiltonian, which, how- 
ever, requires the evaluation of an additional integral with mixed coordinate second 
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derivatives [HI HE]- To avoid this additional integral evaluation we develop here an 
alternative approach. 

In short, our computational strategy in the LFCC formalism to obtain eigenstates 
with various angular momentum quantum number is as follows. In our variational 
procedure we use basis functions, which are eigenfunctions of the total spatial angular 
momentum operators, L 2 and L z [T3]. This is the simplest way to make sure that 
we obtain angular momentum eigenstates, since rotational "contamination" cannot 
be removed by a simple subtraction of a term from the full Hamiltonian [TBIITT]. 
Then, we investigate the effect of the parametrisation of the basis functions on the 
translational contamination of the total energy and correct for it during the evaluation 
of the integrals in the LFCC formalism. 

The paper is organised as follows. Firstly, in Section we present the necessary 
theoretical details. In Section [3j we identify the translational energy contamination 
of the matrix elements. Next, in Section H] we present numerical results for the two 
lowest rotational states of the H 2 molecule. 

2 Theoretical Details 

In this section we present the theoretical details of the variational procedure applied 
to solve the full quantum Hamiltonian without clamping the nuclei. We introduce 
the notation essential to this work. Firstly, we define laboratory-fixed Cartesian coor- 
dinates (LFCCs) and translationally invariant Cartesian coordinates (TICCs) [18] as 
well as the corresponding non-relativistic quantum Hamiltonians. The introduction of 
the TICC formalism is necessary to better understand the mathematical properties of 
the LFCC formalism, which we then use in the calculations. Then, we define the basis 
functions using explicitly correlated Gaussian functions (ECGs) [T9H23] and the global 
vector representation (GVR) [23H26] . Finally, we present the corresponding matrix 
elements and point out the parametrisation problems in the LFCC formalism. 

The generalised eigenvalue problem corresponding to the matrix representation of 
the Hamiltonian is solved using the standard linear algebra library routines of LAPACK 
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(Version 3.2.1) [27] through the Armadillo framework (Version3.4.0) [28] . 

2.1 Coordinates and the Quantum Hamiltonian 

The non-relativistic quantum Hamiltonian for n + 1 particles with m, masses and g, ; 
electric charges is in atomic units 

n+1 _^ n+1 n+1 

?'— 1 i= 1 j>i 

where the vector r = (r l5 r 2 , . . . , r„ +1 ) collects the LFCCs. The translational motion 
of the center of mass has traditionally been separated by introducing some set of 
TICCs [18], x = (xi, x 2 , . . . , x n ) and the Cartesian coordinates of the center of mass 
(CMCC), X CM , by a linear transformation 



x 



TICM 



X 

X 



CM 



[U x ® l 3 )r (2) 



where 1 3 is the 3x3 unit matrix and U x is a non-singular constant matrix with the 
restriction 

n+1 

J2(U x ) tJ = 0, with i = l,2,...,n (3) 

j'=i 

and 

{Ux)*+ij = ™- j /m tat , with j = 1,2, ...,n + l , (4) 

which guarantees the translational invariance for the coordinates x and gives the def- 
inition of the center of mass. The notation m tot = X^"=i m o was introduced. The 
general form of the quantum Hamiltonian corresponding to the new coordinates x and 

x CM HE1E51E2! is 

n n n+1 n+1 

n mi = -i££ - +-a, cm + £ £ , «* , <5) 



i=i j=i i=i j>j iv-'y J / i 



with V^. = (d/dxn.d/dx^.d/dx^) 



n+1 



= ^2( U x)*k(U x ) jk /m k with z, j = 1,2, . . . ,n (6) 



k=l 
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and 



= (U:% - (U- 1 )^ with /../ 1.2 n. (7) 

In Eq. (jSJ) the translational kinetic energy of the center of mass, T CM = — l/2m tot A X(; 
separates, and by subtracting it the translationally invariant (TI) Hamiltonian, if TI = 
H — H CM , is obtained. Mayer and Rokob [30] have recently shown that the same result 
can be obtained by requiring that the total momentum of the system is zero. 

Traditionally the TI form of the Hamiltonian has been used in pre-BO calculations 
[inilUlIIl] • in this work we explore the consequences of the direct usage of the LFCC 
formalism and the original Hamiltonian, Eq. (JT|). 

2.2 Basis Functions 

The basis functions are defined similarly to Ref. [H]. The spatial part is constructed 
using ECGs [T9H23] and the GVR [2UI25]. In what follows we introduce the notation 
relevant for this work. 

The basic structure of a basis function written in LFCCs, r, is 

ij LML (r;A,u,K) = \v\ 2K+l Y LMl (v) exp ("^ T (A <g> l 3 )r) (8) 

with the global vector 

v = (u® l 3 )r . (9) 

The values collected in u = (ui, u 2 , ■ ■ ■ , u n+1 ) are variational parameters. Further 
variational parameters are K e N and the elements of the symmetric matrix A which 
describe correlations between the elements of r. The matrix A must be positive definite 
to ensure that the basis functions are square-integrable and have a non- vanishing norm. 

Furthermore, in Eq. flH]) Y LM (y) is a spherical harmonic function of degree L and 
order M L and v represents the direction of the global vector. Thus, the parity of 
the basis function ip LML is (— 1) L ("natural parity"). The quantum numbers L and 
M L correspond to the total "spatial" (orbital plus rotational) angular momentum 
operators, L 2 and L z . 
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Next, we consider the transformation of the LFCCs, r, to some set of TICCs, x and 
the CMCCs, X CM , in order to better understand the properties of the LFCC formalism, 
which we intend to use in variational calculations. Upon a linear transformation of 
the coordinates the mathematical form of the spatial functions remains the same, and 
only the parameter matrices have to be transformed 

1 



ij LAlL (r;A,u,K) = \v\ 2K+l Y LMl (v) exp [-^r T (A <g> l 3 )r 



\v\ 2K+L Y LMr (v)exp ( --x^ lCM (A 



1 



(x) 



ls^TICM 



with 



V 



(u <g> l 3 )r 
(u {x) ® l 3 )a; 



TICM 



(10) 

(11) 

(12) 
(13) 



According to the definition of the ceticm coordinates, Eq. (j2J), the parameter matrices 
can be transformed back and forth as 



A [x) = UZ'AU: 1 & A = UiA^U, 



(14) 



and 



— (@J), the parameter matrices 



Due to the translational invariance condition, Eqs. 
have block structure [Hj 



and it (a:) 



(s) 



c A 

By combining this property with Eqs. (fT0|) - (flTT) we obtain 



(16) 



tf>LM L (r; A, u, K) = \v\ 2K+l Y LMl (v) exp (~^ T (A ® l 3 )r) (17) 
= Iwl^y^ (0) exp (-i« T (Jl w <8> l s )« - ^X* M ) (18) 



and 



v 



[u <g) l 3 )r 



,0) 



l 3 )cc + c u X 



CM , 



(19) 
(20) 
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which means that the system is at rest only in the case of c A = and c u = [T ^ [3T | [32] . 
At the same time we have to ensure that our earlier requirements for the basis function 
parameters, u and A are fulfilled. While the condition c„ = does not cause any 
problems, we find that Ca = does so. If Ca = the matrix is singular and 
since A and A^ are related by a linear transformation, so is A. This violates the 
requirement for A to be positive definite. Thus, we can conclude that it is not possible 
to parametrise ECGs in the LFCC formalism so that the basis functions are square 
integrable with a non- vanishing norm and at the same time the system is at rest. 

According to the block structure, Eq. ffTBT) and the back-transformation, Eq. fTH|) 
the exponent matrix is constructed in our calculations according to 

/ n+l \ 

= -a l3 (l - <y + V a ik 5 l3 + c A —^--f- . (21) 

\ — / I 'hot I 'hot 

2.3 Matrix Representation of the Hamiltonian 

In spite of the parametrisation difficulties described in Section |2T2"| we intend to use the 
LFCC formalism to construct the matrix representation for the quantum Hamiltonian 
because of its original simplicity. We repeat here only the necessary expressions from 
Ref. [H] and for the original integral derivation see [26] . 

The matrix element of the kinetic energy operator for the 7th and Jth quasi- 
normalised basis functions is Hi 



£,3/4 



1 " <^|V^ (M®1 3 )V^ |V>j 

i=l 

> ?«/ / V 9«j / \ V qu i qu J ) KpujVjPuj.u 




m=0 

Pu I U I , < V ^\ PU J- U J , (7 , o^\ Pu I> U J 



,1 . 



R u + {Kj - m) — ±-L + (Kj - m) J —^- + (L + 2m)- 



HhKjKjm (22) 

PUj,Uj ' PUj,Uj ' PUj.Uj 

where M is a diagonal matrix with Ma = 1/rrij. The H LKlKjm terms are precalculated 
factors 

H 4'»(L + m + l)! 1 m] 

L «i«J™ (K I -m)\(K J -m)\m\(2L + 2m + l)\ y/F KlL F KjL 1 1 
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in order to increase efficiency and ensure numerical stability with the terms 



K 



V- 4"(L + m + l)! 

2-j (K - mV(K - mVm,U?J, + 9,m + 9.V. { ' 



KL 



m=0 



(K - m)\{K - m)!m!(2L + 2m + 2)! 



stemming from the quasi-normalisation 

hM = «V>shM)2 Ze{/,J}. (25) 

Furthermore we introduced short hand notations for terms which we have to study 

in terms of their dependence on c A . One term we have to study is 

= det(2A J )det(2A J ) 

det(A /7 )det(A /J ) 1 ] 

stemming from the origin centered Gaussian functions. Furthermore we have to analyse 
the terms depending on the global vectors 

Pu x ,u Y = u x A "i 1 j u y with X,Ye{I,J} (27) 

Puj.uj = -uJA^AjMAjA^Uj (28) 

Puj.uj = -u^A^AjMAjA^Uj (29) 

P Ul ,Uj = uJA^AjMAjA^Uj, (30) 
the terms from the quasi-normalisation 

qu z = -u^A^Uz with Z e {I, J} (31) 

and the term 

Rjj = ^Tr [AjjAjMAj] (32) 

which can be associated with the radial motion of the system |33j . 

Here, we immediately see that the singularity of A/ and Aj would cause terms in 
Eqs. (|2T)|) -( 13"2"|) to be not defined. The singularity is introduced if the c A = selection 
is made to guarantee translation free expressions. 

Since the potential energy terms in Eq. ([1]) depend only on the inter-particle dis- 
tances, we can choose any c A > value and evaluate the matrix elements without any 
problem and the resulting potential energy matrix elements are independent of the 
value of c A . So, they are not discussed here any longer and are evaluated according to 
Ref. [H]. 
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3 Identification and Elimination Strategy for the 
Translational Contamination 



In Section 12. 2\ we noted that the ECGs take the same mathematical form, indepen- 
dent of whether we choose the LFCCs r, or some a? T icM- This simple transformation 
behavior also transfers to the expressions of the integrals. Due to this property we can 
study the influence of Ca on the different terms in Eqs. f l26|) -f l32|) . only the parameter 
matrices A 7 , Aj and tt/, Uj corresponding to r have to be replaced by their according 
expressions in terms of Aj, Aj and Uj \ Uj corresponding to cc TI0M . The parame- 
ter matrices are related by the transformation given in Eqs. (I14p -(IT5 ]) . and Ay , A^ 
and uf\u^ have block structure, Eq. (TT6l) . 

Firstly, let us consider the R u term of the kinetic energy matrix elements, Eq. (132|) . 
explicitly. We analyze the properties of Ru with respect to c A using the TICC for- 
malism: 



R u = -Tr [AjjAjMAj] 



-Tr [{A u ) )-'AfU x MU T x Af ) ] . 



and also exploit the block structure of the matrices: 



(33) 



(x) 



c A 



with z E {I, J, IJ} and V ' X MU\ 



(x) 











(34) 



The block structure of U X MU^ follows from the conditions of translation invariance, 
Eqs. dm) - ®, for x. We can thus separate the c A dependent terms in R u : 



Ru = Ru + \caCm = ^Tr [(^-^V^] + \caCm 



The c M factor of the linear contribution is 

c M = (U x MUl) n+hn+1 
and thus with c M = 1/ m tot 



n+l 



(35) 



(36) 



R 



ij 



R]1 + l C - 



4 m 



(37) 



tot 
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Next, we investigate the c A dependence of the remaining terms and factors, Eqs. (126]) - 
(|3ip . of the matrix representation of the kinetic energy in the LFCC formalism, 
Eq. (122]) . In Eq. ( 126]) we find that any contribution of c A cancels. In Eqs. (|3~T]) and ( 127]) 
we see that c A only contributes if c u > 0. Since we require c u = 0, this contribution 
is eliminated. Finally, in Eqs. (l2"g]) - ([3"0"]) we find that the contribution of c A to the 
exponent matrices cancel. In short, only the R ItJ term, Eq. (132]) . has a non- vanishing 
(linear) c A dependence in the kinetic energy matrix element, T u . 

Thus, if c A = was chosen, the translational dependence vanishes, but the expo- 
nents matrices, Aj, Aj, are singular without having an inverse. Thus, for an imple- 
mentation in a computer program we can choose a non- zero value for c A , and eliminate 
the translational contribution explicitly by subtracting 3c A /(4m tot ) from the R u ma- 
trix element. This is a simple computational strategy which we are going to follow in 
the LFCC formalism. 

We also note here that the direct variational optimisation of all parameters, includ- 
ing c A here, would be another option that has been suggested already in the literature, 
e.g., in Ref. [2]. From the theoretical details presented so far, we understand that the 
total energy with c A > is always an upper bound to the total energy free of the 
overall translation of the system, E tot (c A ) > E Ti . But as we have shown for c A = 
several parameter matrices incorporated in E tot (c A ) are singular, so in spite of the fact 
that the limit exists, the application in a computer program is problematic (c„ = is 
chosen throughout the discussion). 

This explains our preference for the approach developed here, which releases the 
translation-free condition for the basis functions and corrects for the translational 
contamination in the kinetic energy explicitly for each basis function, I = 1, 2, . . . , N. 
The details of our algorithm are: 

1. Generate, optimise or read in the a Iti j values for i = 1, 2, . . . , n + 1, j = i + 1, i + 
2,...,n + l. 

2. Construct the elements of the exponent matrix in the LFCC formalism as 



10 



/ n+1 \ 

(Aj)^ = n,.„(l - <y + > a J)ift <5 y - + c A 

\ / ''''tot '''tot 

with i,j = 1, 2, . . . , n + 1 and > 0. 

3. Due to the c A > choice the matrix Aj is non-singular, and 1/ det(Aj) and 
A7 1 can be calculated. At the same time the total kinetic energy contains some 
translational contamination. 

4. The translational contamination is eliminated by replacing R u , Eq. (132)) . with 
Ru — 3c A /(4m tot ) in the expression of the kinetic energy matrix element, T u , 
Eq. d23). 

Throughout this computational strategy for the elimination of the translational con- 
tamination in the LFCC formalism we have c u = 0, c A > 0. 



4 Numerical Examples 

In this section we present numerical applications using the LFCC formalism. The 
appearance of the translational contamination and its elimination according to the 
strategy described in Section [3] are demonstrated. 

Our test cases are the lowest energy levels of the para-H 2 (L = 0,p = +1) and the 
ortho-H 2 (L — l,p — —1) molecules both in the singlet electronic state. These are the 
two lowest-energy rotational states of the hydrogen molecule. L is the total spatial 
(orbital and rotational) quantum number and p is the parity. The wave functions 
are obtained by a direct solution of the linear variational problem using 1500 basis 
functions with an optimised parametrisation taken from [2]. Here we use the LFCC 
formalism exclusively, thus all parameters were transformed first to the LFCC repre- 
sentation, according to Eqs. fn^ - ffT5]) . The exponent matrix for each basis function 
was constructed according to Eq. fT2T]) . During this transformation we were free to 
choose c A to simulate different levels of translational contamination, while c u = was 
fixed. We used the same c A value for each basis function. Table [I] collects the results 
of our numerical calculations. 
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[Table 1 about here.] 

The first column of Table [I] lists the five values of c A which have been studied. 
The next two columns provide the total energies of the system with and without 
translational contamination, respectively. The last two columns list the kinetic energy 
contribution with and without translational contamination, respectively. We observe 
that those energies, which contain translational contributions (columns 2 and 4) de- 
pend on the value of Ca and that an increase of the value of Ca causes an increase of 
the energy, according to Section [3j We also note that the corrected energies (columns 
3 and 5) are independent of c A . Furthermore, we list the translational correction given 
in Eq. flBTI) . Hence, these results give a numerical confirmation that we have identified 
and eliminated the translational contamination depending on the c A basis function 
parameter, while choosing c u — 0. 

5 Conclusion 

In this work, we used for the first time laboratory-fixed Cartesian coordinates (LFCCs) 
with an explicit correction strategy for the translational contamination in pre-Born- 
Oppenheimer variational calculations with various angular momentum quantum num- 
bers. Our work was motivated by the inherent simplicity of the LFCCs and the corre- 
sponding operators. 

Instead of transforming the coordinates we accounted for the translational and 
rotational invariances of the isolated many-particle problem by using an appropriate 
form and parametrisation of the basis functions in the variational procedure. The basis 
functions were constructed using explicitly correlated Gaussians and the global vector 
representation, and as an extension of our earlier work [2] we focused here on the 
usage of LFCCs and the problem of translational invariance. 

First of all, we observed that it is impossible to parametrise explicitly correlated 
Gaussian functions (ECGs) in such a way that the total system is at rest in LFCCs 
and at the same time the basis functions are square-integrable with a non-vanishing 
norm. 
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Fortunately, it was possible to devise a simple computational strategy to circum- 
vent this problem. So, for the sake of the stability of the numerical computations we 
released the translational constraint on the basis functions, by choosing a non-zero 
value for the free parameter, ca, in the LFCC parametrisation of the ECG exponents, 
and then explicitly corrected for the translational contamination in the kinetic en- 
ergy integral expressions. This correction term is a simple constant, which depends 
linearly on Ca- Its form was derived by considering a few mathematical relationship 
of the formalism: a) the properties of the linear transformation between LFCCs and 
translationally invariant and center-of-mass Cartesian coordinates (TICMCCs); b) the 
corresponding transformation of the basis function parameter matrices; c) the fact that 
the parameter matrices are block diagonal in TICMCCs. 

It was also shown that the uncorrected total energy E tot (c A ) is an upper bound 
to the translation-free total (intrinsic) energy. Thus, in principle, we could obtain 
this value by the variational optimisation of the LFCC parametrisation (implicitly 
including the c A value in the optimisation). We prefer, however, our explicit treatment 
for the elimination of the translational contamination, because in the c A = limit the 
exponent matrix, A, of each ECG would be singular, and thus the numerical evaluation 
of 1/ det(A) and A" 1 would be impossible. 

Finally, to demonstrate the numerical applicability of our approach we calculated 
the lowest two rotational levels of the singlet hydrogen molecule corresponding to the 
para and ortho proton spin states, respectively. 

The presented LFCC formalism with the explicit translational contamination cor- 
rection is an alternative but equivalent to the traditional approaches using some set 
of translationally invariant Cartesian coordinates with the Cartesian coordinates of 
the center of mass explicitly separated already in the Hamiltonian, e.g. [T0j[T3]. The 
simplicity of the LFCCs is an appealing choice for the variational solution of the 
Schrodinger equation with the non-relativistic Hamiltonian. Furthermore, one can 
think of more complicated operators for which the usage of the simplest possible co- 
ordinate representation and the avoidance of any coordinate transformation is more 
than just a comfortable option. 
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Table I: The appearance and elimination of the translational contamination of the total pre-Born- 
Oppenheimer energy using the laboratory-fixed-Cartesian-coordinate (LFCC) formalism. The ground- 
state energies, in Eh, of the singlet, para-H2 (L = 0,p = +1) and ortho-H2 (L = l,p = —1) molecules 
are calculated using the basis function parameter set of Ref. |14] , The corresponding translationally 
invariant energies are E TI (L = 0) = -1.164025 026 E h and E TI (L = 1) = 1.163485 167 E h , respectively. 61 
Results corrected for the translational contamination of the LFCC formalism are indicated with the 
"corr." subscript. 







^LF.corr. 


(^ LF |T LF |^ LF ) C 


(* LF |T LF |* LF ) c c orr . 




para-H 2 L = 


0.01 


-1.164022985 


-1.164025026 


1.164027045 


1.164025004 


0.000002041 


0.10 


-1.164004614 


-1.164025026 


1.164045416 


1.164025004 


0.000020412 


0.50 


-1.163922966 


-1.164025026 


1.164127064 


1.164025004 


0.000102060 


1.00 


-1.163820906 


-1.164025026 


1.164229124 


1.164025004 


0.000204120 


2.00 


-1.163616786 


-1.164025026 


1.164433245 


1.164025004 


0.000408240 



ortho-H 2 L = 1 



0.01 


-1.163483125 


-1.1634* 


55167 


1.163487203 


1.1634* 


55161 


0.000002041 


0.10 


-1.163464755 


-1.1634* 


55167 


1.163505573 


1.1634* 


55161 


0.000020412 


0.50 


-1.163383107 


-1.1634* 


55167 


1.163587222 


1.1634* 


55161 


0.000102060 


1.00 


-1.163281047 


-1.1634* 


55167 


1.163689282 


1.1634* 


55161 


0.000204120 


2.00 


-1.163076927 


-1.1634* 


55167 


1.163893402 


1.1634* 


55161 


0.000408240 



a The parametrisation of the internal basis set (V) was taken from Ref. |14j . The number of basis 

functions was 1 500, the order of the polynomial prefactors was 2K £ [0, 20]. For the proton-electron 

mass ratio the value m p /m e = 1836.152 67247 was used, and thus mtot = 3 674.305 34494 m e . 

b Free parameter of the explicitly correlated Gaussian basis functions expressed in LFCCs. The same ca 

value was used for each basis function in the basis set. The value c u = was used throughout the 

calculations. 

c Elf , ^lf : Eigenvalue and eigenfunction of the full Hamiltonian expressed in LFCC using the basis 
function parameters (V,ca)- 

Tlf : The kinetic energy operator expressed in LFCCs. 

"corr." : Correction for the translational contamination in the LFCC formalism, as explained in 
Section [2j 

<^Tr. = 3cA/(4mtot) : Translational contamination of the kinetic energy, see Eq. ()37p . The corrected, 
translation- free value is obtained as i?LF,corr. = -Elf — ^Tr. • 
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